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ABSTRACT 


The objective of the present work is to predict the 
Vibrational behaviour of a coupled structure Cformed by 
coupling two identical beam elements), when the response 
characteristics of each individual beam member is known 
from experimental modal testing. The beam, elements are 
connected at their ends using a single bolt. This arrange- 
ment facilitates the two beams to be bolted together at any 
desired angle between the beams. The out-of-plane oscilla- 
tions are considered so that when the longitudinal axes of 
the beam members are not collinear then the bending and tor- 
sional modes get coupled. Theoretical and experimental 
analysis is carried out. Theoretical analysis gives a good 
prediction of the natural frequencies of the coupled structure 
for various angles of coupling. Impedance coupling method is 
used in the experimental analysis. The substructure response 
data are obtained using single degree of freedom circle fit 
method. The coupling models considered account for the 
coupling of beams at various number of points with various 
degrees of freedom Ctranslational and rotational). The 
experimental models have failed to yield a good prediction of 
the response of the coupled structure. 



CHAPTER 1 


INTRODUCTION 


1. 1 Introduction : 

The dynamic testing is being extensively incorporated for 
the design development, quality control and qualification of 
products in Industries. Experimental ^fodal Analysis is an 
important technique which is being widely used for this purpos* 
Much development has taken place in the last two decades in 
this area of applied science. The primary objective is to 
determine or predict the dynamic behaviour of the test object. 
This information can then be used for fault diagnosis, design 
improvement and evaluation of operational capability of a syste 
Modal testing is also used to pinpoint and Improve the modellir 
inaccuracies of finite element analysis. 

An important application of experimental modal analysis 
is in the prediction of response of a coupled structure. The 
coupled structure is an assembly of substructures. If the 
response characteristics of these substructures in isolation 
is known, then the response of the coupled structure can be 
obtained using a coupling technique. The prediction of respons 
of a coupled structure from the responses of substructure offer 
the following advantage. The substructures can be modified in 
isolation and the combined effect of these modifications can be 
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studied quickly without having to perforn the modal testing 
and analysis for the complete coupled structure. Some typi- 
cal applications in the aircraft structures are listed below. 

Ci) "todifi-cation of mass, stiffness and damping proper- 
ties of aircraft com.ponents Cii) Modification of Inertia 
properties in the usually large number of various wing-v;ith- 
stores configuration of fighter aircrafts. Ciii) Attachments 
of additional boosters to a launcher system. Civ) Payload 
alterations of space systems, for instance, by exchanging the 
entire payload module of a shuttle. 


1.2 Literature Review : 

This section restricts the literature review to coupled 
structure responses using modal testing. The literature conce 
ing the development of Modal testing is covered in recent 
works by Jaspal Singh [1] and Vasudeva Rao [2]. 

Coupling techniques can be broadly classified into impe- 
dance coupling method and modal coupling method [3], Modal 
testing combined with finite element models is employed widely 
for analysing the coupled structures [4]. 

In impedance coupling method, the measured frequency 
response data of substructures of a coupled structure are used 
to predict the response of the coupled structure. The number 
of degrees of freedom that are to be considered at each 
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coupling point is an important aspect. This depends upon 
the symmetry of structure and nature of the coupling. The 
coord5.nates associated with degrees of freedom can be classi- 
fied as translational coordinates and the rotational coor- 
dinates. "^easurem.ent of oscillations along translational 
coordinates is straightforward and can be done, for example, 
v;ith an accelerometer. The techniques available for the m.eas- 
urem.ent of rotational coordinates and the techniques for pro- 
viding moment and torque excitation are limited. 

Smith [ 5 ] has developed a method for the application of 
couple. Twin shakers are used to apply either a direct force 
(with the two shakers in phase ) or a couple (with them in ant 
phase), Ewins and Sainsbury IS~1 have outlined a method which 
uses a single shaker. A specially designed exciting block is 
attached to the test specimen at the coupling point. Forces 
are applied and measurements are made at different points and 
in different directions to obtain the mobility matrix which 
considers the rotational coordinates. 

Sattinger [73 discussed a method in which rotational 
mobilities are derived as the spatial derivatives of the 
translational mobilities using finite difference technique. 
Chen and Cherng [83 have given a simple method for the measure 
ment of rotational mobility. Rotational effects are approxi- 
mated by reducing parallel force systems separated by a small 
distance to an equivalent force-couple system^. 


K 


The literature available in the area of coupling with 
modal test data is relatively snail. Ewins and Sainsbury [6] 
have analysed the coupling of a beam, with a block. Impedance 
coupling method was used. Two components were bolted together 
at a single point. Three degrees of freedom were considered 
at the coupling point. The results exhibited a good prediction. 

An interesting problem, of analysis of coupled helicopter 
carrier has been carried out by Ewins, Silva and .Maleci [9], 
Goyder IlO] has discussed structural modelling from measured 
frequency response data. An error analysis of impedance coupllr 
method is made. It is concluded that the response predicted frc 
Impedance coupling method can only be best treated as a guide tc 
the average response of the structure rather than as a detailed 
prediction. 

Ewins C3ll describes two modal coupling methods. The modal 
models of the substructures are used to determine the required 
substructure impedances. In method 1, the impedance coupling 
technique is recast in a form which directly gives the solution 
for modes of the coupled system. Method 2 (flexible coupling 
method} assumes the existence of a spring and/or dashpot element 
between each pair of connected subsystems. Both the methods 
finally reduce to a standard eigenvalue problem, 

A discussion of different types of modal coupling namely 
Rigid, Flexible and Mixed coupling can be found in [113* Convex 
gence problems due to frequency range truncation and Incomplete 
load conditioning in the coupling points, statically determinate 
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and overdeterminate coupling conditions are also discussed. 

’'•lost of the work found in the literature employ a combined 
modal testing and finite element technique to analyse the 
coupled structures. A good review of literature of such a 
combined technique Is given by Grief Chen and Cherng 183 

have outlined a combined technique. Coupling of a beam with a 
beam and a beam with a plate was discussed. Lee [123 has used 
modal testing, condensation method Cof PEM) and sensitivity 
analysis for dynamic modelling of structures with bolted and 
bearing j oints . Condensation technique [static and dynamic 
condensation) is a numerical technique used for reducing the 
computation by reducing the order of the matrices [133. 

Component Mode Synthesis is another reduction technique 
used in finite element modelling of the components. The basic 
principles involved in component mode synthesis is explained 
in [1^3. Grief [43 has given a review of the literature in 
this area. 

1.3 Objective and Scope of the Present Work : 

This work was undertaken to Investigate the coupling of 
two beam elements. Two beam elements bolted together by a 
single bolt form the coupled structure. This type of coupling 
facilitates the beams to be coupled at different Included 
angles. The objective of the work is to examine the variation 
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in the vibrational behaviour of the coupled structure when 
the two beams are coupled at different angles. The objective 
was to develop both theoretical and experimental formulations. 
Both theoretical and experimental analysis are carried out under 
free-free conditions. 

A closed form theoretical formulation is presented in 
Chapter 2. This formulation considers the coupling of out-of 
plane bending vibrations with the torsional vibrations of the 
coupled structure when the beams are coupled nonsymmetrically . 

Impedance coupling technique is used to study the coupling 
experimentally. The number of points considered for coupling 
and the associated coordinates at each point are explained in 
Chapter 3 - 

In Chapter 4, the modal testing is used to obtain the 
modal parameters for the substructures and these are used to 
predict the response of the coupled structure. A comparison 
is also made therein to compare the predicted and experimental! 
measured responses. Chapter 5 presents a general discussion 
on the results obtained from various coupling models. 



CHAPTER 2 


THEORETICAL FOPJWLATION 


2 . 1 Introduction : 

This chapter deals with the theoretical formulation 
for the coupled bending and torsional vibrations of the coupled 
structure (section 2.2). In section 2.3, a comparison of theore- 
tically obtained natural frequencies with the experimentally meas- 
ured ones is given. 

In the theoretical formulation, the coupled out-of plane 
bending and torsional vibrations is dealt with. The theory of 
bending vibrations is based on the Bernoulli-Euler beam theory. 
Literature dealing with the coupled out-of plane bending and 
torsional vibrations is not known. Chang [173 discussed the in-pZ 
vibrations of frames and trusses with inclined members. 

2.2 Formulation for Coupling of Out of Plane Bending and 
Torsional Modes of Free Beams 

Consider two free-free beami elements of rectangular 
cross-section joined at their ends (Fig. 2.1). Figure 2.1 
shows the global CX-Y - Z ) coordinate system and the local 
CXi - - z^} coordinate systems attached to the beams. The 

included angle between the two beams is denoted by 2r. Let 
yj^Ct, x^) and ^j^Ct, Xj^} be the bending and torsional response 






q 


coordinates. The equation of motion for and 9^ are. 



+ B‘ 


ib 


y. 

d 


= 0 


and 


d 




d^ 0, 


d x2 


0 . 


where. 


62 = ~i . g2 

ib Pj. ^ ®iQ 



^i ^pl 


2.1(a) 


2.1(b) 


and is the modulus of elasticity, is the area moment of 
Inertia, A. is the area of cross-section, p. is mass density, 
is shear modulus, is torsional constant [16] and is 
the polar moment of inertia of the ith (i = 1,2) bar. (Geometric 
and material properties used in the present work* are given in 
Appendix A) . 

To study the natural modes, the oscillations are 
assumed to be harmonic so that 


Ji = y^Cx^lSlnpt 

§1 Ct,x^) = 9^U^)Slnpt , ^^ 2 , 


C2.2) 


Substitution of (2.2) in (2.1) yields the form of solution as 
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YlCxi) = d^SinK^^x^+d2CosK^l^x^+d-Sinh K^^x^+d|^CoshK^^x^ 

(2.3) 

Y 2 CX 2 ) = d^SinK2^X2+d^CosK2^X2+d^SinhK2|^X2+dj^CoshK2^X2 


= ‘^9®-" ’^19 + “*10 %e ’'i 

= d-,1 Sin K2Q X2 + d., 2 K2Q X2 . 

in whj-ch.. 


C2.4) 


SQRT ( — ) and (i=l,2) 


p is the frequency of vibration in radians/sec. The constants 
of integration, d^s are to be determined from the boundary 
conditions . 


The bending moment and torsion at cross-section at 
distance x^ are. 




d X? 
1 


(2.5) 


and 


T. 


= G^C. 


dXj^ 


Ci = 1,2) 


( 2 . 6 ) 


Free'- free coupled system shown in Pig. 2.1, has the 


following six natural boundary conditions. 
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At = 0 


do^ 

d^ 


d^ y- 


d x: 


= 0 


d’ y- 

d Xt 


= 0 


At X 2 = I 2 


dft, 

C 

dx^ 


d^y. 


dx; 


= 0 


= 0 


d^y, 

dXn 


= 0 


(2.7) 


Following are the six compatibility conditions at the 
joint (x^ = 1^ and X 2 = 0): 


(i) Displacement of both the beams should be equal 

Yl = 72 


( 2 . 8 ) 


Cii) Shear force balance at the joint yields 


% ^1 


dx 


Yi d-^y^ 


dX; 


(2.9) 


Ciii) Rotation of beam 1 and beam 2 In X-dlrection (Global) 
should be equal (Fig. 2.2) 
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= © 2 ^°) Cos’? 




n¥ 


dx. 


or 


, , dyp(o) 

©lCl]_) = -© 2 '^°) c:os 2r - — Sin2r (2.10) 


(iv) Rotation of beam 1 and beam 2 in Z -direction (i.e, 

Global) should be equal (Figure 2.2) 


d y^Ci^) 

g— 


dypCo) 

© 2 ( 0 ) Sin'i' + — — Cos¥ 


or 


dy^^Cl^) 

33r 


© 2 C 0 ) Sin 2T 


dy 2 (o) 

— ^ Cos2r 


( 2 . 11 ) 


(v) Moments in Z-direction on the beams should balance 

(Pig. 2.3). 


or 


Ml = T 2 Sin'? + M 2 Cosf 


d^y . (1 ) d95(o) 

®1^1 " ^2^2 ~dxj 


d^yp(o) 

E_I- 5 — Cos2r 

^ dx^ 

( 2 . 12 ) 


(vi) 


or 


Moments in X-dlrection on the beams should balance 
(Fig. 2.3). 

Tl = T2Cos'i' - M2 SinT 

■30,(10 a9,(0) a^y,(o) 

TxT 1x7 - hh -ri— ^^"2? 

1 2 dXg 


(2.13) 



*1 = *1 
X2= 0 



FIG. 2*2 RESOULITION OF COORDINATES IN x AND z DIRECTION 


*1 'M 

X2 



z 


FIG. 2-3 RESOLUTION OF COUPLES IN x AND z DIRECTION 



The first boundary condition of (2.7) j l.e. 9|(o) = 0, 
gives one of the constants of integration dQ = 0. Substitution 
of C2.S) and (2.4) in (2.7) to (2.13) gives a set of homogeneous 
equations of the form 

IM] {d} = {0} (2.14) 

where 

[M] is llx 11 square matrix and {d} is column matrix 
of eleven constants, such that 

{d} - { d^ jd2 jd^ ,djj ,d^ ,dg ,dy ,dg ,d^Q ,d|j^^,d^2 ) 

and the elements of matrix [M] are given in Appendix P. 

For non-trivlal solution, the deteriment of matrix 
[M] should vanish. The parameters and in the elements 
of matrix [M] contain the frequency, p. The value of p 
at which the determinant bectsmes zero is noted as the 
natural frequency of the system. 


2.3 Verification of Theoretical Results : 

Two beams, of dimensions and properties outlined in 
Appendix A were coupled at various values of the included 
angle. Equation (2.14) was used to determine the natural 
frequencies of the coupled structure. Numerical computation 
was carried out at the step of IHz. 
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An experimental verification for these results vjas 
also performed. Two beams of the same specifications were 
bolted. An impact hammer was used to excite 

‘ . coupled beam structure under free-free condition. 

Table 2.1 gives a comparison between natural frequencies 
obtained theoretically and experimentally for two values of 
included angle C2r). The torsional and bending modes for 
2r = 180 ° are identified in Table 2.1 since these modes are 

uncoupled. For other values of included angle it will be more 
meaningful to talk in terms of symmetric and unsymmetrlc modes. 
Included angle being zero is a special case. Pig. 2.^‘ presents 
the results for r = 0 but these results are not clubbed with the 
other results obtained for different values of 2r, because the phys. 
interpretation can not be attempted without looking at the modes. 
Fig. 2.5 Ca) and (b) give two experimentally measured frequency 
response functions of the coupled structure for 2r = 90° and 
2r = 180° respectively. 

In the experiment, the two beams overlap on each other 
reducing the effective length of the coupled structure (in 
case of 2r = 180°1- But in the theoretical formulation the 

beams couple at one point making the length of the coupled struct- 
ure Cin case of 2r = l80°) equal to the sum of the lengths 

of Individual elements. Moreover, the nature of coupling between 
bending and torsional modes is strongly dependent on the included 
angle. Therefore, it cannot be said conclusively that the natural 
frequencies, as obtained from experiment, will always be higher 
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Table 2.1 Theoretically and Experimentally Measured Values of 
Matural Frequencies for 2r = 90 and 2r = 180° 


Included 
angle (2r) 

Theoretically obtained 
natural frequencies 
(upto 2000 Hz) 

Experimentally meas' 
natural freauencie 
(upto 2006 Hz) 

o 

o 

286 

272 


310 

304 


532 

608 


544 

624 


848 

964 


1032 

1004 


1537 

1612 


1563 

1624 


1879 

— 


180° 

89 (bending mode) 

100 


246 (bending mode) 

288 


48l (bending mode) 

528 


6I8 (torsional mode) 

704 


795 (bending mode) 

908 


1187 (bending mode) 

1196 


1230 (torsional mode) 

1316 


1(^57 (bending mode ) 

1808 


1845 (torsional mode) 



Frequency (<*>) 
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than the theoretical ones. But generally the experimental values 
are on the higher side. 



CHAPTER 3 


THEORY OP IMPEDANCE COUPLINC USING EXPERIMENTAL PRPs 


3 . 1 Introduction : 

This chapter explains briefly the theory of impedance 
coupling methods. Different types of Frequency Response 
Functions (PRPs) and the relevant definitions are detailed 
in Appendix B. First a method is outlined for the coupling 
of two substructures at a single coordinate. Then the pro- 
cedure is outlined for coupling at more than one coordinate 
with more than one degree of freedom. Inclusion of the 
coordinates in the analysis which are not involved in the 
coupling is also examined. Section 3.3 explains in detail the 
various coupling models used in the present work. 

3 . 2 Impedance Coupling Method : 

This method is also called as Dynamic Stiffness Method [33. 
The underlying principle is given below. 

Consider two components A and B (Pig. 3.1) which are 
coupled at a single coordinate to form the coupled structure C. 
The coupling coordinate is designated by x^ and the coincident 
coupling points on substructures A and B are denoted by numeric 

1 « 






UJ 


o 

cr. 



O 

Ll 
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The coupling point on the coupled structure ’C', is also 
denoted by 1. 

If the system A at the connection point is excited harmo- 
nically by e^“^, the resulting response x^ is given by. 


^a ^iojt _ a / % f.a lot 


or 


a: 




similarly for subsystem B, 


X- 


= a (cj) . f. 
^ 1^1 


(3.1) 


(3.2) 


(3.3) 


Applying the compatibility and equilibrium conditions at 
the connection coordinates one gets. 


X 


a 

1 




(3.4) 


f 


c 

1 


.a 


C3.5) 


Substitution of C3.2) and (3.3) in ( 3 . 5 ) and using the 
equation (3,4), the relationship between the receptances is 
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or 




.a 

XiXi 






= Z' 




^ 1^1 


( 3 . 6 ) 


( 3 . 7 ) 


where Z denotes the impedance. 

Thus the PRF o? the coupled structure *C* can be obtained 
in terms of the FRF properties of the substructures analysed 
separately. 

The above analysis can be extended to the systems connected 
at more than one coordinate. At each coordinate more than one 
degree of freedom may be considered. The coordinates which 
^re not involved in coupling may also be included. The whole 
process involves the partition of the subsystem impedance mat- 
rices into submatrices which are involved in coupling and those 
which are not. The submatrices Involved in coupling are added 
to yield the resultant Impedance matrix. Therefore, if Co^U)] 
and [Z Cw ) 3 are the receptance and impedance matrices, respect- 
ively, of subsystems A and [a^Cto)] and [Z^Cw)3 are the corres- 
ponding matrices of subsystem B, and if the coupling coordinates 
are collectively designated by q, the coordinates not Involved in 
coupling by p for subsystem A, then the receptance matrix is 
partitioned as 


24 


[a^Co))] 


X z 

p p 


x z 
q p 


a 

z 

D q 


z z 
q qi 


C3.8) 


IZ^U)] = [a^Cto)]"^ 




z^ z^ 

X X 

q p q q 


(3.9) 


Similarly, th.e impedance matrix for the subsystem B is 
given by, 

[z^(a)) ] = Ea^(a))]"^ = 

where ’'r* denotes the coordinates on the subsystem *B' which 
are not involved in coupling. 

The order of each submatrix in C3.9) and C3-10) depends 
upon the number of coordinates involved in the coupling and 
the number of degrees of freedom considered at each coordinate. 

Now as before, applying the compatibility and equilibrium 


X X 

r r 


X X 

q r 


X X. 


X X 


q Qj 


(3.10) 


condition, one can obtain. 
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1 

Ea^Co))! 


© 


1 

ta^Coj)! 


( 3 . 11 ) 


where, © sign indicates the addition of partitioned matrices 
Involved in coupling. 

In oth-er words, the impedance matrix [Z (o))] for the 
coupled structure can be obtained in the form 


rz^'Ca))] = e rz^(co) ] 


C3.12) 


X X 

P P 


,a 

■'x X 

q. p 


X X 
r r 




X X 

P q 




Z®' +z^ 

x„x„ x„x„ 

q q q q 


(3.13) 


'J'h.e recept^nqe matrix Ea^Cco)J for the coupled system is 
obtained by the inverse of EZC(u)3 as 


Ca^Ca)) 3 = IZ^'Ca))) 


-1 


(3.14) 


Prom modal testing point of view, the receptances required 
for th.e substructures may be evaluated by any one of the options 
mentioned below. 
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CD Using the modal Model and evaluating the receptance matrix 
using the formula 

[aCoj)] = [0] r (Ap - ^ (3.15) 

where 

xl = 0 )^ (1 + 

C2) Calculating the receptance FRF matrix obtained by direct 
measurement . 

C3) Calculating the regenerated receptance FRF matrix from 
the measured FRFs. 

The results presented in Chapter 4 are by using the last 
method of utilising the regenerated FRFs. 


3.3 Number of Points and Degrees of Freedom Considered for 

Coupling : 

The equilibrium equations for coupling of substructures, 
which are undergoing bending oscillations, will require the 
conditions on shear forces and bending moments. Therefore, 
two degrees of freedom, namely, transverse displacement and 
angular rotation should be assigned to the coupling point. 

The simplified coupling analysis can start with coupling of only 
translational degree of freedom. 
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Figure 3.2, shows the type of coupling being considered 
in this work. The two beams are coupled with a bolt and nut. 

The simplest coupling model is assumed to couple the two 
beams at only two points. This is later increased to coupling 
at four points. 

Thus the different combinations of number of points 
selected and the degrees of freedom assigned to each point is 
expected to yield differing predictions for the coupled structure 
The various models cosidered in the present work are presented 
in the following sections. 

The notations used for the coordinates and points are as 
indicated in the corresponding figures. Initially the coupling 
is explained for an included angle of 2r = l80°. 


3 . 3.1 One-point, One-dof Coupling ; 

Figure 3-3 shows the coupling points and the degree of 
freedom considered on each substructure A and B and also on the 
coupled structure ’C’. Prom eqn. (3*6), the receptance of 
coupled structure will be. 


_1 ^ 1 

“yiYi “yiYi 


b 

‘^ 1^1 


If the substructure A and B are identical, then the substructure 
cL Id 

receptances a„ ,, and „ will also be similar. This will 

Vl ^1^1 

lead to the coupled structure receptance as 
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_1 ^ 1 ^ 1 

a'' 

YlYl 


where. 


a 

‘YlYl 




Yiy 


i-^i 


1 a 

2 


or 


YlYl 


= i 

2 y^Yi 


which means that the coupled structure will have a similar FRF 

at the coupling point as its substructures. In fact the response 

of the coupling point on the coupled structure and the response 

of the corresponding point on the substructure will show peaks 

at same frequencies. For other natural frequencies of the 

coupled structure, this coupling point happens to be a node and 

hence will not show any additional peaks. Obviously this 

^1^1 

model is applicable only to situations like longitudinal oscilla- 
tions of free beams. 

3-* 3. 2 Two-point One-dof Coupling : 

JPigure 3.4 shows the coordinates and the degrees of freedom 
considered for the substructure and for the coupled structure. 
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Prom equation (3,11), 


Ca""] 


© 


[ct^] [ab] 


i2>^^ ^ “yo>y 


( 3 . 16 ) 


where , 


[a^3 


y2>y3 



a 

“^2^2 

a 

“y2y3 


CO 

on 

b 1 

a 


a 

^"^2^3 

a 

”^3^3 


b 

“y3y2 

b 

a 

^2^2 

tan 


where. Maxwell's reciprocity theorem has been Used to make 
the receptance matrices as symmetric matrices. The impedance 
matrices are obtained as, 


CZ^.] 


^2^3 



^3 


[Z^] 


^3^2 


= [a^] 


-1 

y3,y2 


Now while coupling the substructures it can be noted 
from Figure 3.4, that the point 2a couples with the point 3 b . 
This means that the impedance matrix of coupled structure will 
be. 



''3b ''2 b 



F16. 3-4 TWO-POINT ONE-dof COUPLING MODEL 
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z^, + z?, , 

^2^2 ^3^3 

^2^3 ^2^3 

Z®’ + 

Z®- + z^ 

y2y3 ysys 

^3^3 ^2^2 


(3.17) 


Tfie two coordinates considered does not explicitly account 
for the effect of rotation due to bending but indirectly the 
effect of slope is introduced. 


3 . 3.3 Four-point One~dof Coupling : 

Figure 3.5, shows the coupling coordinates considered 
for the substructures and for the coupled structure. Procee- 
ding in the manner explained in section 3 . 2 , one can arrive 
at. 


rz^i 



,yg,y3,y5 



^ 3^3 


Symmetric 



a 

n-a 

^3^2 

z" 3 

^3^6 

a 

^5^5 

z ^ 

^5^2 

Z^ 


Z®" 

^ 2^2 

Z®" 



Z^ 

.b 

'^25^6 

Z^ 

z"^ 

y2y5 

,b 

'^ 6^6 


Z^ 


^3^3 

Z^ 

^3^5 

>> 


Symmetric 
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Noting the fact that the points 335 * 2,6 of beam A get coupled 
vjith the points 2 , 6,33 5 of beam B, the impedance matrix of the 
coupled structure can be written as. 


[Z*"] 


y3,y5,y2»y6 


yd. +7^ ’7^ +Z^ 

‘^y3y3 ^y2y2 y3y5 

yd. +7^ 

Symmetric 



Z^ 

+Z^ 




a ^2^ 

z^ 

+Z^ 

^ 5^2 ^ 6^3 

^ 5^6 ^ 6 ^ 

+z^. 

Z^ 

+Z^_ _ 

^2^2 ^ 3^3 

^ 2^6 ^ 3 ^ 


( 3 . 18 ) 

This model can account for coupling between bending and torsional 
oscillations . 


3.3.4 One-point Two-dof Coupling ; 

While considering th.e bending oscillations, the eqili- 
brium conditions in coupling the two beams can be written in 
more precise terms if rotational coordinate is also assigned 
to the coupling point. This procedure is explained below. 

The coupling coordinates considered for the substructures 
and for the coupled structure are shown in the Fig. 3.6Ca). 

The rotational receptances required are calculated by reducing 
a linear parallel force system applied at two close points 
into a single force and moment system £81!. Fig. 3.6Cb) shows 
the principle involved. 







3 ? 




1 


«* *1 ♦ »2 
M*= f,l4*l 


FIG. 3-6 (b) PRINCIPLE INVOLVED IN INTRODUCING ROTATIONAL 
COORDINATE 
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The receptances associated with rotations at the coupling 


point CH are point receotances, namely, a„ „ and a ^ and 

^l^i ^i^i 

cross-point receptance, namely a„ q . These are obtained in 


i'^i 


terms of substructure point receptances Ca„ ^ and a 


, ) and 

transfer receptance Ca The receptances a etc. can 

^ 1^2 ^ 1^2 

be measured directly using a single accelerometer. For this 

reason the receptances associated with rotatational dof i.e. j 

a„ c etc. are derived in terms of measured translational recep- I 
^ 1^1 ■ ! 
tances a etc. as explained below. [For definition of various ' 

^ 1^2 ' 
receptances [ 3 ] see Appendix BJ . I 


If the distance ax CFig. 3.6Cb)) between the parallel lines j 
of action of forces f 2 _ and is small enough, the system of 
forces approximately reduced to an equivalent 

force-couple system at either point. If the equivalent force- 
couple system is considered at point 2, then the force at point 
2 will be, 


f* = 

C3.19Ca)) 

and the couple will be 


M* = f]_CAx) 

C3.19Cb}) 


The rotation, at point 2 can be expressed as. 


0, 


■ ^1 - 
Ax 


C3.20) 



The 


Thus the two coordinates at point 2 are and Q^. 

objective Is to obtain th.e receptances a and a so 

^ 2^2 ^ 2^2 

that the equivalent force and the equivalent moment can be 
used to write, 


^2 " ^ “^ 2^2 ^ 


( 3 . 21 ) 


^2 = “ y 2©2 " "© 2^2 


The eqn. (3.19) can be used to express and In terms of 

* * 
f and M as 


\ = M /AX 


( 3 . 22 ) 


f2 = f 


Ax 


on the other hand and y 2 can also be written as 




C3.23) 


0 „ .. + 0, f, 

2^2 


^2 “y^y^ ‘1 ' "y,y, ‘2 


Then eqns. G-201, G-. 22) and G.231 can be arranged In the 
form of eqn. _G.21) such that, 
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^2 


a 

^ 2^2 

Ot ot 

^ 1^'2 ^ 2^2 

A X 

1 

c\j 

a* 


a -a 

^ 1^2 ^ 2^2 

a,, ~2a 

YlYl y2^2 ^1^2 


Ax 

— 1 

OJ 

X 

<1 


(3.24) 


or {y> = Cal„ ^ {f} 

«y 2 > ^2 


(3.25) 


Figure (3. 6 (a)) shows the coupling of two substructures 
at points C2a, 3a.) and C3b, 2b). This is considered as 
equivalent to one point coupling at point 2a and point 3b 
where two .coordinates (i.e. translation and rotation) are 
associated with these points. In order to keep the same dire- 
ction for rotation at point 3b and point 2a, we define, 


©. 


b 


y2b - y3b 


and 0, 




and M 








( 3 . 26 ) 


where the subscripts 'a* and ’b’ refer to the substructures 
A and B respectively. 

Eqn. C 3 . 26 ). are used to express the receptance matrices 
for substructures A and B in the form of (3.25) as 



hi 


* 




... 


'* ^ *1 

i'a 

.®a. 

a 

= [ci^] 

.”aj 

and 

^b 

- [o'"] 

* 

L t) _ 


where , 


[a^] = 


¥2^2 


8. Bi 

”^2^3 " “^2^2 


a 


(3.27) 


fY 8 8 

^ - a 


^2^3 


^2^2 


n,®’ J- ™9- o 3- 

®TT ,T ■*■ “,T ,T “ 2a„ 

^ 2^2 ^ 3^3 ^ 2^3 


"a 


(3.28) 


[a*>.] 


y^yj 


b b 

“y2y3 “ “y3y3 


b b 

«„ - 0i„ 

^2^3 ^3^3 


c'= + »»> - 2a‘’ 

^2^2 S'3y3 7 yjyj 


The impedance matrix of the coupled structure (Pig. 3.6(a)) can 
he obtained in the form (from equation 3.12) 


[ 5 =] 


ys.Sj 




^8 


+ z 

='3='3 

+ 

^2^2 

jy U 

y3e 


Q + 

y2«2 

"3® 


(3.29) 



As an Illustration for this model, the tip receptances. 


for the undamped free-free beam vibrations given by Bishop 
and Johnson [I 5 ] are used. The example cited here is for tne 
coupling of tvfo identical beam elements of properties detailed 
in Appendix A. The theoretical receptances and the receptance 
mat-rix for- the coupled structure are given in Appendix C. 

Figures 3.7(a), (b), (c) show the substructure (Appendix 

‘ X, Fig. C.l) accelerances a,, , a„ ^ and a^ ^ respectively. 

^1^1 ^l‘=^l 

The material daiuping is assumed to be zero. Figures 3 . 8 (a) and (b) 

show the accelerance af and a2 ^ at the coupling point of 

^ 1^1 ^ 1^1 

the coupled sti'ucture (Appendix C, rig. C.2), respectively. These 
figures display accelerances rather than receptances. This is to 
maintain the uniformity with the results presented in Chapter 4, 
where the conversion from receptance to accelerance is also expl- 
ained. 


Another example is taken where the material damping is 
also considered. The formulas for calculating the point and 


transfer receptances (i.e. a , a and « ), as given 

*^ 2^2 ^ 2*^3 ^ 3^3 

by Snowdon [18], are given in the Appendix D. Coupling of two 


Identical beams of Appendix A is considered and the rotational 


receptances (consisting of a_ ^ )and <* 00 ) are ootained in the 

^ 2^2 ^ 2^2 


form as given by equation C3.27). Then equation (3.29) is used 
to obtain the receptance for coupled structure. Figs. 3.9(a), (b) 


and (c) show a , a and a of the substructure accelerance 

^2^2 ^ 2^3 ^3j3 c -c 

and Figs. 3. 10 (a) and (b) show (^y y , ag g ) and ay g respective' 

2 2 00 00 




2^2 







1^1 


Pig. 3.7(c) * 
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FIG. 3-tl {q) 






FIG. 3-11 (b) TWO-POINT TWO-dof COUPLING MODEL 
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'’•.■.' 0 - 10 '" Coupling; : 


"^.llCa) and (b) shov/ the coupling coordinates 
■"“.nr ‘ di'V-ij 'cn bhe substructure and the coupled structure. The 
para 'lie: force systers acting at points (la, 2a), C3a, 4a) on 
?ul r.'zr^ucTur® A, (Ih, 2b), (3h, 4b) on substructure B are reduced 
■; r '■■r. ■‘•ci' valent system of force and couple at points la, 3a on 
A .an I '^r and 4b on B respectively. Figure 3*llCa) also shows 
the expression 'or the equivalent force and moment. With these 
relati'^ns, the ^’eceptances involving rotation are written in 
terr.r, of receptances ^or translation only. The procedure is as 
fn soction ^.3.4 and the obtained receptances are. 
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yi,y, ysyii 


a 
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(3.30) 


Other elements can be obtained using the symmetricity of the 
recentance matrix. 


Elements of [a*^] 
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/ifter this the impedance matrices of the 
calculated. These are coupled to get [Z°Cw)]y^^Q 


substructure are 

^ as 
i.y3.93 


CZ°U)] 




Cz^(u>)]y^^Q^^y3^e3 0 [z ^“5Jy2,©2*y4,ei^ 

C3.'32) 


And finally the receptance matrix for the coupled structure is 
obtained by the inverting the impedance matrix of this coupled 
structure from (3 •321. 
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FIG. 3-12 (a) 
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FIG. 3-12 (b) 








~ r.t '^^h'^^ee-do'^" Coupling : 

1' tors’’ onal vibrations vjere also found to be present 
■ r. the structure, then the equilibrium conditions in coupling 
the t'A’o beams can be written in more precise terms if rotation 
associated with torsional vibration is also assigned to the 
coupling point. This procedure is explained below. 

Figure 3.12(a) and (b) show the coordinates considered 
■^or reducing the parallel force system acting at points (la, 2a, 
3a) and (lb, 2b, 3b) to an equivalent system of force, moment 
and torque at points la and lb on the substructures A and B 
respectively. Figure 3* 12(c) show the coupled structure and 
the coordinates associated with the coupling point. 

Using the relations of equivalent force, moment and 
torque as given in Fig. 3 •12(a) and (b), the substructure recep- 

tance matrices and ta^(a))]y^ 3 _,ei^,©l,j can 

be obtained in a similar way as explained in section 3.3.^. This 
yields. 


Elements of Ca^I]y-j_,93_^jjj©]_i];i: 
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(3.34) 


These are used in the impedance coupling method for obtaining 
the receptance of the coupled structure as given below. 
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CHAPTER a 


EX^ERP-^EHTAL RESULTS AND DISCUSSION 

'4 . 1 Introduction: 

In this chapter, a brief explanation of the experimental 
setup and a procedure of testing adopted are given in section 4.2. 
The section 4.3 presents the predicted response of the coupled 
structure. Each subsection of 4.3 gives the results obtained from 
different coupling models. Substructures are analysed separately 
and the measured frequency data obtained from substructures are 
used to predict the response of the coupled structure. This is 
checked against the experimentally obtained FRF for the same 
coupled structure. 

The FRP^s displayed in this chapter are for accelerances 
Ca) but the modal parameters extraction are done from the recep- 
tance Ca). The accelerance is measured instead of the required 
receptances. Receptances are calculated using, 

A 

2 

Accelerance, aj^Coi) = - - ui • 

This additional computation required is due to the fact that 
the double integration process in charge amplifier results in 
the usage of high gain in amplifier which Introduces noise. 

This can be clearly seen by observing coherence for both 
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'icclerar.ee and receptance measurement [1]. 

'*’he measurements m.ade on the test specimen Csubstructure 
or coupled structure) is for the free-free condition. The test 
specimen is hung over a frame using light strings to simulate 
the free-free condition. 

^ ^ Experimental Setup and Details of Testing : 

The scheme of experimental setup is shown by the block 
diagram in Pig. 4.1. The instruments used in the present work 
are described briefly in the folloxAring. 

Impact excitation using a hammer is a commonly used 
technique. Structures of masses 2 kg to 1000 kg can be excited 
by impact hammer. An impact hammer consists of different types 
of tips and heads for obtaining different frequency and force 
level ranges. Below the tip is a force transducer for detecting 
the magnitude of force of excitation. An impact hammer type 8202 
CBriiel and KjaerX is used in the present work. 

Piesoelectric miniature accelerometer CB&K, type 4374) 
is used to pickup the response. An accelerometer is preferred 
over other types of pickups because of its smaller size and 
wider frequency range of operation. The miniature pickup has 
a mass of only 0.65 gm and with small dimensions (base diameter = 
5 mm). It is suitable for measuring signals at closely placed 


locations . 



HIOKI Recorder 
8801 


Block Diagram of Experimental Setup 
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Two charge amplifiers, one for response signal (B&K, 
t,/pe 2<^^5) and the other for excitation signal (B&K, type 2626), 
qre used to amplify weak signals into signals strong enough to 
be measured by the analyser. 

The spectrum analyser receives two time signals, 
namely, the excitation and the response signals. The analyser 
operates on these signals to determine the system transfer- 
function or frequency response function. The operation performed 
on the tv/o signals depends upon the function selected on the 
panel of the analyser. Analyser of type SM - 2100B (Iwatsu) is 
used in the present experimentation. Coherence function is 
selected as the operating function. The use of coherence function 
automatically calculates the auto-power, cross-power, coherence 
function and the transfer functions. These are registered in 
different memory blocks by the inbuilt functions. 

The selection of frequency range and the data length 
decides the frequency resolution. Once the frequency range is 
selected, the data length can be selected depending upon the 
frequency resolution required. 

In the post- processing (while extracting modal parame- 
ters) of the PRPs, the PRP data registered in the analyser is 
transferred to the IBM-PC/XT via GPIB (IEEE-4887 bus. The data 
can also be sent back from PC to analyser. Por taking hard copy 
print -outs of PRPs displayed on the CRT of the analyser, a 
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'■e”:ory Hi corder Type B801 (Hioki) thermal plotter is used. 

Some of the precautions to be taken while fetching the 
olcnnls nre noted below. The strength of the signal should be 
V. *thln the scale selected for the A/D converter of the analyser, 
"'■nee the strength of the signals are set using amplifier, the 
;../D scale can then be set properly so that the signals lie 
within this scale. If this is not ensured then the digital data 
stored in the analyser memory will be faulty, i.e. the data gets 
clipped at the limiting value of the A/D converter scale. 

A flat spectra of the force signal upto the frequency 
range of Interest should be ensured. This can be achieved by 
proper selection of the tip of the impact hammer, by maintaining 
a correct duration of impact and by avoiding multiple hits. Plat 
nature of the autopower of the force signal indicates a correct 

impact . 

The noise content of the signals should be kept minimum. 
By monitoring the coherence function, which should have a value 
of unity for good quality of signals, this can be achieved. High 
noise content in the region of ahtiresonance can be neglected 
as the data in this region is of not much importance in most 
of the experimentation. 

Some of the limitations of the analyser may be noted. 

The unavailability of zoom facility causes difficulty in the 
extraction of moda parameters of lightly damped structures or 
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0 ' "repkly excited nodes. If multichannel facility were 
:iv^-lnole, tnen a good amount of measurement time can be 
saved especially when measurements at large number of co- 
ordinates are needed as in the case of two-dof and three- 
do coupling models. The Ixvatsu analyser provides only 
two types of windows namely rectangular and banning windows. 
Rectangular window is used for both excitation (impulse) 
signal and response (transient) signal. Exponential, window [ 3 ] ' 
xfould have been better for guiding the response signal as there 
is no windowing effect when a transient signal Is guided through 
a rectangular window. The Iwatsu analyser does not permit the 
selection of independent windows for two channels. 

^*3 Experimental Response of the Coupled Structures Vs . 

Predicted jResponse~TIsing Impedance Coupling wHerT^he 
Substructures are Identicar i 

The coupling models discussed in Chapter 3 have been 
applied to the coupling of two identical beam elements. The 
substructures and the coupled structure along with their 
dimensions are shown in Pig. 4 . 2 (a) and (b). The two beams 
are of MS and coupled at their one end using a bolt. A small 
washer was Introduced between the two components to reduce the 
contact area within the experimental limits. 

The response predicted using the various coupling 
models is compared with the directly measured experimental 
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resrr'nse fron the coupled structure. The results were first 
cr-r.p"red for an included angle between the two beams of 180° 
(?y = l3o°). At this angle the torsional and bending vib- 

rations are uncoupled. 


l^.^.l Steps Involved in obtaining the Predicted Response of 
Coupled Structure : 

The first step is to extract the modal parameters 
from the experimentally obtained PRPs of substructures. These 
modal parameters are used to regenerate the receptance of the 
substructures. The regenerated receptances are then used for 
predicting the response of coupled structure. Since the two 
substructure beams are identical, measurement on one beam is 
sufficient as long as the coupling points on both beams are 
on identical locations. The associated steps are better illus- 
trated through examples. Thus the steps involved in obtaining 
the regenerated receptances ^ s (corresponding to trans- 
lational coordinates only) are explained in section 4. 3 .2.1 
and the steps involved in obtaining the regenerated a^^g^s 

1 J 

(which corresponds to combined translational and rotational 
coordinates) are explained in section 4. 3 *2. 3* 
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^ . 2 Predicted Response Through Different Coupling; Models : 

Tn the following sections the method involved in 
rredicting the response of coupled structure using various 
models is explained. 


. P . 2 . 1 Two-Point One-dof Coupling : 

With reference to Pig. 4.3, the coupling points on 
substructures A and B for this model are assumed to be the 
points 2 and 3. The coupling is assumed to be at (A 2 , B^) 
and CA^, B 2 ). The receptances to be used in equation (3 . lH are 
the regenerated receptances. The procedure involved in obtain 
ing the regenerated receptances is explained in the following. 


The first step is to measure the substructure recep- 
whlch 

, /comprises of the point and transfer recep- 


tance Ca^Cw)] 


72 , 


tances measured at point 2 and 3- The single degree of freedom 
circle fit method II] is used to extract the modal parameters 
from the experimentally obtained receptances. These modal 
parameters u)p, rij,) are then used in the equation 


a 


N 

Z 

r=l 



r^Jk 

27 

- 0) + in CD } 


C4.1) 


to regenerate the receptances. The regenerated receptances 
are used rather than the measured ones to avoid the noise 
present in the experimentally measured receptances. 
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^Igure shows a measured and regenerated 


y^y.^ 


s-ibstructure P^^s for a frequency range of 1 KHz. Figure 
■'.5(a) and ^ig. 4.5(b) show a„ „ and a„ „ measured and 


y2y:^ 




regenerated substructure PRPs for a frequency range of 
5 KHz. It can be seen from these figures that the regen- 
eration is satisfactory and the modal parameters extracted 
are acceptable. 


After obtaining the regenerated receptances, the 
next step is the calculation of substructure Impedance 
matrices [Z®'(a))T. and CZ^(a))]„ „ as a function of 

frequency using equations (3.9) and ( 3 . 10 ). With the sub- 
structure being identical, the regenerated receptances 
of one beam are used for the other beam also. Therefore with 
the coupling points 2 and 3 In Pig. 4.3, the receptance matri- 
ces for beams A and B are related as: 


If Ca^] 



then [a^] 



(4.2) 


The complex inversion involved in obtaining Impedance from 
receptance (or receptance from Impedance) is carried out 


using the method given in Appendix E. The Inversion is obtained 


at each frequency step over the frequency range considered. 


The substructure impedance matrices are then coupled 

according to the equation (3.17) to obtain the Impedance matrix 

[Z‘^Gj>)] of the coupled structure. The receptenace matrix 

^ 2*^3 
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V coupled is then obtained by inverting 

' 2>''3 

Using the procedure outlined above, the predicted 
response was obtained for frequency range of 1 KHz. Fig. 
H.f(a) shows the predicted accelerance and the Pig. 4. 6(b) 
shov;s the experimentally measured response of the coupled 
structure for 1 KHz frequency range. It is clear that the 
results of the predicted response are totally unsatisfactory. 
Moreover, the predicted receptance of coupled structure shows 
the peaks at the natural frequencies of the substructure 
(see Pig. ^1.4). Measurement of response of the substructure 
over a wider frequency range results in inclusion of more 
number of modes and is expected to Improve the results. 

Pig. 4.7Ca) shows the predicted accelerance when the measure- 
ments were made upto 5 KHz Cwhich includes upto 8 modes). It 
is observed that there is no Improvement in the predicted 
response and the peaks still occur at the substructure natural 
frequencies. Pig. 4.7Cb) shows the measured accelerance for 
the coupled structure for 5 KHz range. 

H.3.2.2- Four" Point One-dof Coupling : 

With reference to Figure 4.3, the coupling points 
on substructure A for this model are assumed to be the points 
2,6,3 and 5. The corresponding coupling coordinates on sub- 
structure B are assumed to be the points 3,5,2 and 6. The 




;nurl‘.nt" is assumed to be at (A 2 , (Ag, B^), (A^, B 2 ) 

:nd (A^, 

The regenerated receptances [a (co)] 


and 


are obtained in the similar way as expl- 
ained in section ^i. 3.2.1. After this the impedance matrices 
C^Coj)]^ ^ ^ and CZ^Ca))]„ ^ are calculated. 

These impedance matrices are coupled according to equation 
(3.18) to get ^ ^ ^ , which is then Inverted to 

<y 2 5 ^ ^ 3 ^ 

obtain [a^Cto)]^ ^ ^ • 

y2»^6»^3»^5* 

Figure 4.8 shows the measured and regenerated acce- 

lerances a . for a frequency range of 2 KHz. The modal 

parameter values of measured „ are listed in Table 4,1* 

Figure 4.9(a) shows the measured and Fig. 4.9Cb) shows 

the predicted a^ of the coupled structure. A careful 
3^2^ 2 

study of Pis* ^•9(b) sliov/s small blunt peaks at around 250 Hz^ 
450 Hz and 1200 Hz which happen to be the 2^^^, 3^*^ and 6^^ 
modes of the coupled structure. But since all the substructure 
natural frequencies have reappeared in the predicted response, 
the prediction is still not a good representation of the 
coupled structure response. 


4. 3. 2. 3 One- Point Two-dof Coupling; 

In this model the coupling is assumed to be on a 
single point (A^, B^) CPlg. 4.3) and two coordinates ®2 

are associated with this point. Initially the receptances 



Table 4.1 : Modal Parameter of Substructure (Torsional Mode Considered 
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V ^ \7 

I Q 3 3 


are measured at 2 and 3. These are 


3 } Ct^ 

.hen T’e venerated as explained in section 4. 3. 2.1. The 


I'esenerated PRFs are used to obtain Ca^Cui)],. ^ and 

j 2 ^ ^2 


r-b 


a (a))1r Q explained In section 3.3.^ of Chapter 3. 

j 3 s 

The distance between the parallel lines of action of 

forces is taken as 25 imi- a and [oi^Cw)3„ q 

^2’^2 ^3’ 3 

are inverted to give [Z^'Cco)] a and tZ Cw)]^ a • These 
impedance matrices are then coupled according to equation 
(3.29) to get [Z°Ca))l„ a which Is then Inverted to get 
[a'°(w)] • Table 4.2 shows one set of sample calculation 

y 2^5 ©2. 

at 0) = 356 Hz (first natural frequency of the substructure). 


Figure 4.10 shows the measured (torsional 

mode,a) = 1242 Hz, is not considered In regeneration) and 

Fig. 4.11 show predicted a^ ^ and q for 2 KHz frequency 
range. Eventhough a direct comparison between 
Cexperlmentally measured) and “9^9^ (predicted) is not 
true Indication of quant-itative differences, the Inferences 
about natural frequencies can be made. From the Figure 4.11, 
it can be noted that there is no improvement in the predicted 
response, as all the substructure natural frequencies have 
reappeared very prominently. A very small peak can be 
noted near the second natural frequency of the coupled 


structure. 


Iitle : ■'’alues of Receptances and Impedances of 

Substructures and Coupled Structure at 
a Frequency of 356 Hz 

(The distance between two points 2 and 3 
= 0.025 m ) 


Ci) Regenerated receptance ; ^ 


PRP 

Real part 

Imaginary part 


1.654x10"^ 

1.905x10"^ 

“^2^2 
a tr TT 

1.915x10"^ 

1.925x10"^ 


5.155x10”^ 

3.815x10"^ 

^3^3 



n 2» 

Sj and CZ®C»)]y^^B 

2 

Cal Receptance 


FRF 

Real part 

Imaginary part 


1.654x10"^ 

1.905x10"^ 

^2^2 

PZ 

_7 

ot ^ 

1. 046xl0”^ 

8.000x10 ' 

^2®2 



Q 

^2 2 

4.765x10"^ 

2.994x10 ^ 


Continued, 
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(iii) 


■'•.7 (Continued) : 


(b) J”:pedance 



^eal part 

Imaginary part 

’^ y , y2 

2.583x10^ 

■1.638x10^ 

-3.019x105 

7 . 221 x 102 

02^2 

1 . 501 x 1 o 2 

-9.587x10 

[“”<^“>1 fl 

and 

• 

0 

!’^3 

(a) Receptance 


FRP 

Real part 

Imaginary part 

a 

^3^3 

“^3®3 

9 

^3^3 

5.155x10“^ 

1. 296 x 10 ”^ 
4 . 765 x 10 “^ 

3 . 816x10”^ 

7.5^5x10"^ 

2.994x10"^ 

Cb} Impedance 

FRP 

Real part 

Imaginary part 

^^3^3 

^383 

^®3®3 

2 . 583 x 10 ^ 

- 6 . 621 x 10 ^ 

3.197x10^ 

- 3 . 019 x 105 

8 . 272 x 1 o 5 

-3.207x10^ 


Continued, . , , 
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( i V ) 


(to)] 




and ra'^Cto)] 


^2 


,©2 


(a) Impedance 


FRR 

Real part 

Imaginary part 

2 .y -- 

5.166x10^ 

-S.039x105 

^2^2 




-6.785x1o3 

8.994x103 

r\ 


4.698x10^ 

-4.166x10'^ 


(b) Receptance 


FRF 

Real part 

Imaginary part 

a 

1.186x10"^ 

1. 165x10"^ 



2.498x10"^ 

1.528x10“^ 

^2^2 

“©2®2 

1.689x10'“^ 

1.254x10"^ 
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7'v-^o-Polnt Two-dof Coupling: 


In this model the coupling is assumed to be at 
and CA^, B 2 ) (Pig. 4 . 3 ) and the coordinates (y]_j9]_) 
and (y^, 9^) are associated with these points. The distance 
between the parallel lines of action of forces (required 


^or calculating tfie equivalent moment} is taken as 10 mm. 


cL 

rz (a))L^ p. „ Q and [Z“(a))I]., q v. q obtained by 

y]_»»xsy3 >^3^ y2,>^2»^4»®4 

inverting ra^(^i:]„ « v ft Cci^(u}]„ ^ v, ft, » 

yi3yijy3s»3 y2»®2»^4»®4 

respectively. These impedance matrices are then coupled to 

get the Impedance matrix of the coupled structure which is 

inverted to get the receptance matrix of the coupled structure 

according to equation C3.321, 


Some of the measured and regenerated substructure 

PRPs are superimposed in Fig. 4.12, for a frequency range 

of 2 KHz. The effect of torsional mode on the receptance 

is very small as this is very weakly excited for this set of 

points 1,2,3 and 4 lying on the center line. Therefore, 

the torsional mode (natural frequency = 1242 Hz) is not 

considered while regeneration. This can be seen by looking 

at the regenerated PRPs in Pig. 4.12, which does not show a 

peak at 10 = 1242 Hz. Table 4.3 gives the modal parameter 

values for the PRPs shown in Pig. 4.12. Pig. 4.13 Cal and 

4. 13 Chi shows substructure a^ ^ and a^ ^ 2 KHz frequency 

^ 2^2 - ■- ^ 2^2 

range as calculated from equation C3.24}. The general nature 
of these PRPs is still same as those shown in Pig, 4.12. 



Table 4,3 ‘ Modal Parameters of Substructure 
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‘'’Igure (b) and (c) show predicted af, q » aS 

_c ^ 3^1 ^3 

and ^ of the coupled structure, respectively. A 

3 3 

snail peak at around 250 Hz corresponding to 2 mode of 
coupled structure has appeared In the predicted response. 
But the prediction Is not good as all the substructure 
natural frequencies have reappeared very prominently. 




Pig* 4.4 t ^^ 373 ^ 

Pig* 4.4 * Measured and regenerated substructure PRP 
(1 KHz range) 



Pig* 4*5* Measured and regenerated substructure PRP (5 KHz range) 



.L^— i:,— i. 








4- f |, « , X...-W . I ii l |T— — i-j 


Fig*4*6(a}t |a^ { 


^2^2 


Pig. 4. 6(b) t laf V I 


72^2 


Fig. 4 .6* Two-point One-dof coupling* Predicted and measured 
response of coupled structujre (1 KHz range) 


A^jTDOQ ■•E0'.-. 




4.7(a) » |a° 1 

yy 


Fig. 4.7(b)* ja^l 


Pig. 4 .7* Two-point One-dof coupling* Predicted and measured 
response of coupled structure (5 KHz range) 







Fig. 
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Fig. 4.10* |a^ y 1 


ED'^:0O0-JDB- ■ . r-aO’* tlOQ- DB- 

Fig. 4.11* |a^^y^|#la|^©^l 


J.4.10 & Pig. 4,11 * One-point Two-dof cc^^ling* Meas^d ^ 
^ ^ predicted response of cotipled strixiture 
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Fig. 4.12 * 1 
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Pig. 4.12 I 

Pig. 4. 12 I Keaaxired and regenerated substructure PRPs 


t! 


Pig.4.13{a) . Flg.4.13(b). 

Rotational ^PRPs calcxalated from regenerated substructure 
PRPs 




Fig. 4.13 * 
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ria, 4.14 Cc) t li| Q { 

3 3 


Predicted responses of 



CHAPTER 5 


GENERAL DISCUSSION 


The original aim of this work was to predict the 
response of an assembly of identical simple beams when 

i) The receptance properties of the simple beam is 
experimentally available, and 
ir) The out of plane bending oscillations are under 

study where the coupling between bending and tor- 
sional oscillations can take place due to unsymme- 
trlc configurations of coupled structures. 

The simplest coupled structure was formed with two 
identical beams with included angle equal to 180°. The 
bending and torsional modes are uncoupled in this configur- 
ation. The attempt to predict the response of this coupled 
structure failed since the predicted response exhibited 
peaks at all the natural frequencies of the substructure 
sJjnple beam. Inclusion of more number of modes in the 
receptance data did not improve the results. Then two point 
coupling with two degrees of freedom associated with each 
point was tried. This improved the results marginally. Lack 
of time prevented trying out more number of coupling points. 
The experimental data then becomes so large that it is really 
not very practical to use this method of prediction. 
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”i, cures ii,12 and Pigs. 3*9 s.re reproduced here in 
n.l to observe that the experimental PRPs match quite 
veil v.’.'th the theoretical PRP (using Appendix D) . As 
evident from the theoretical results in Figs. S-lOj the 
impedance coupling method employed is perfectly correct. 

Then the anomaly could only result from the fact that the 
jiature cf impedance coupling is such that a small Inaccuracy 
in the receptances of substructure causes Inaccurate results. 

In order to narrow down the source of error, a trial 

was made with coupling of two nonidentical beams. The beams 

vrere similar in all dimensions and properties except that 

the length, of one beam was Increased to 0.42 m. Fig. 5 *2 (a) 

shows the two substructures and the points considered for 

coupling. The "one point two degree of freedom coupling” 

(Section 3.3.4) was applied, the coupling point being the 

point number 1. Pig. 5.2(b) shows the regenerated accelerances 

a®’ and a^ . The predicted accelerance for coupled 
^ 1^1 

structure is shown in Pig. 5.2(c) which also includes the 
accelerance obtained experimentally for the coupled structure. 
The predicted response is still not good but it is evident that 
the peaks in predicted response do not always correspond to the 
substructure natural frequencies. It is felt that if the 
difference in lengths of the two beams was larger, the predicted 

closer to the expected response. 


response would appear 












^,y. 


E^qperimental 


Pig* 5.1(b) t Theoretical and Eaperlaental Phase Plots for 

COtop arisen 





Pig. 5.2(b) 

Pig. 5.2(b) » Regenerated s\ab structure PRPs of beam A(1 * 0,3in) and 
beam B (1 « 0.42 m) 


Predicted. 

Pig.5.2(c)* Cfee-point Two-dof 
coupling* Predicted and measuired 
coupled structure PRP. 




Pig. 5.2(c) 
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An error analysis of the Impedance coupling of two 
identical beams is expected to reveal the exact source of 
error in predicting the response of coupled structure. 
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APPENDIX A 


GEOMETRIC AND MATERIAL PROPERTIES OP BEAM ELEMENTS 


Beams considered are of mild steel material. These 
have the following properties: 


E 

V 

G 

P 

a 

b 

1 

A 

I 

I 

] 

C 


Young's modulus = 2.1x10^^ N/m^ 

Poisson's ratio = 0.3 

Shear modulus == 2 (l+v ) — ~ 8.1x10^^ N/m^ 


Mass density = 7.8x10^ kg/m3 

Width of beam = 0.05m 

Thickness of beam = O.OO 6 m 

Length of the beam = 0.3 0.42 m depending 

on the case 

ii p 

Area of cross section = 3x10 m 


3JO 

Area moment of inertia = 


3 


12 


Polar moment of inertia = 
Torsion constant 


ab' 


k.a^b^ 


12 


9x10“^° m^ 

a^b 

12 


(a^+b^) 


where, k for rectangular cross section depends on the a/b ratio [I 6 
for a/b = 8, k = 0.312. 

In the present case a/b ^ 

=> C = 3.32x10“^ 
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APPENDIX B 

FREQUENCY RESPONSE FUNCTIONS (FRPs ) 

Frequency response function (transfer function) Is the 
ratio of output of a system t^o Its Input. Depending upon the 
output selected, the FRFs are classified Into three groups 
as explained below. 


(1) Accelerance [5(ii))(] : Accelerance Is the ratio of acceler- 
ance to the force of excitation. 


aCio ) 


acceleration Cx) 
force (f) 


Cll) Mobility [YCoi)] : Mobility is the ratio of velocity to 

force of excitation 


YU) 


Velocity (x) 
force (f). 


(ill) Receptance taU)] : Receptance Is the ratio of displace- 
ment to force of excitation 


UUl] 


displacement (x) 
force 


For harmonic vibration, the relationship among the above 


mentioned FRFs is given below: 
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aCuj) = = -( 13 ^ a (w) 

* 

= p- = io) a(a)) 

For a beam, the responses and X 2 at two points 1 
and 2 In terms of receptances can be obtained as 

Xi = f3_ + a^2 ^2 

^2 ^ «12 “22 ^2 



The various terminology used with the FRFs is explained 

below: 

A point receptance is one where the response coordinate 
and the excitation coordinate are identical. A transfer 
receptance is one where the response and excitation coordinates 


are different . 



9 ^ 


These are further subdivided into direct and cross- 
receptances. Direct receptance is one where the type of 
coordinates for response and excitation are Identical - for 
example, both in y-direction. p^oss recentance is one where 
the type of coordinates for response and excitation are not 
rdentical like one ^s y-^direction translation and other is a 
- oment in the z-direction. Following are some typical examples 

Ci) Direct point receptances^ 

_ Response in y-dlrectlon at 1 
“y^yi Force in y-direction at point i * 

_ Rotation about z-axis at point 1 
“©i^i ~ ' Moment about z^-'axis at point i 

Cii3 Cross point receptances, 

a = Translation in y-direction at point i 
°'^i®'i Moment about x or z-directlon at point 1 

Ciii) Direct transfer receptance, 

_ Translation in y-dlrectlon at point i 
Force in y-direction at point J 

Civ) Cross transfer receptance, 

_ Translation in y-dlrection at point 1 
“y^ 9. Moment in x or z-direction at point J 


Qh 


14 


Cli' 

“ A p 





El 


A 


Area 

of cross section. 

m^ 

P 

= 

Mass 

density, kg|m 3 



= 

Frequency in rad/sec 




SinXil 

Slnhx il 


^3 

=r 

(CosXil 

CoshXJl )-l 



= 

CosXil 

SinhXil - SinXJl 

CoshXJl 


Fg = CosXil SlniUJl + SlnXX. CostiXH . 

Tiie receptance matrix for the coupled structure is given by 


or. 


Ca°Ccoi J 




+ Pf) 




2F5Pg . 


El X. F, 


-Pj-/X^ 0 

5 


For the maximum value of or ^ , the denomlnatoi 

F.Fg + F2 ^1^1 

of the factor tj ~ should be minimum. Substituting 

2r g • JcLl A • 2 

for F^, Fg, P5 and Fg and simplifying. 




+ P?) 


CCosXJl. CoshXJl + 1) 


2Fg Fg. . F^ 


EIX.CCos 2 xs,.Cosh, 2 X]i - ll 
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APPENDIX D 

ONE POINT 2 DOF COUPLING USING THEORETICAL RECEPTANCES FOR A 

DAMPED S TRUCTURE 

The point receptance at point ’0* and the transfer recep- 
tance points p and 0 as given hy Il83 are. 

Cl + li) Cn* a5 e* 


(1 + ji) (n a5 
“op = - ,2 

where, 

w - frequency of excitation 

0* = CCSh.OCSh.C)^ + CCh.Cl CCh.C)^ 

- CSh.SKSh.S),, - CCh.S)CCh.,S),, - 1] * 

^ Cn al 

IC5h.C+l)C5h.C-Ch.Sl + C Ch.C+1) (Sh.C-Ch.S)3 * 

^ ^ Cn al 

ICe*+6*l Cosh n*Xa + C e* -B ) Cosn*Xa 

*- C 0 *+ K*+ -x*)Sinh n*Xa + Ccf* - k* + x*)Sinn*XaJ 

in which. 



V 
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APPENDIX E 

IN^rERSION OF A COMPLEX MATPTY 

If the inverse of a complex matrix [[A] + i [b]J is given 

by [[C] + i[D]3, then, 

([A] + i[B]) dC] + 1 [DJI = t 1 


( [A] EC] - 

CB] 

CD] 

1 + i (IBJ CCJ + IA3CDJ1 = c 1 : 

] + iIo3 






CE.l) 

Prom 

CE.l) 




• 

[A] 

\ 

1 — 1 
o 

! i 

1 — 1 

L-J 

t — 1 

= L 1 3 

M.zi 

CB] 

CC] + 

f — i 

> 
t — r 

CD] 

= CO] 

CE.3) 


Solving (E.2) and CE.3) for [C] and [D] we get, 

[C] = C CA] + CB] [A]"^ [B]}"^, 

[D] = -CA]"^ [B] C [A] + CB] [A]“l [B])-l. 

The real matrix inversion required above, was done using 
usual inversion process by Gauss-Jordan's method. 
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APPENDIX P 


Elements of Matrix 


M(l,l) 

= 

MCI, 3) = 0 

MCI, 2) 


-®1 Kib 

MCI, 4) 

= 


MCI, 5) 


MCI, 6) = MCI, 7) = MCI 

MC2,1) 

= 

- Eili 

MC2,2) 

= 

0 

MC2,3) 

= 


MC2,4} 

= 

MC2,5) = MC2,6) = MC2 
MC2,11) = 0 

MC3,1) 


“®l^l ^Ib 

MC3,2l 

=: 

®1% ^Ib ^^’^^Ib^l 

MC3,3) 



MC3,4) 

== 

E^I^ Sinh l-]_ 

MC3,5) 

= 

V2 ^ib 

MC3,6) 

= 

0 

MC3,7) 

= 

-^2h 4b 

MC3,8) 

=: 

MC3,9) = MC3,10) = 

MC4,ll 


-Eili SinK^^l;^ 

MC4,2> 

= 


MC4,31 


Sinh 

MC4,4) 


®1^1 ^Ib ^Ib ^1 


= M(l,10) = M(l,ll) 


= MC2,9) = MC2,10)= 
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0 



mC4,6) 

= 

"®2^2 ^2b 



M(4,7) 

= 

0 



mC4,8) 

=: 

Vz 4b 



mC4,9) 

=: 

0 



MC4,10) 

= 

-®2 ^2 '^29 Sin2r 



M(4,ll) 

= 

0 



MC5,l) 

- 

M(5,2) = MC5,3) = M(5,4) 

= M(5,5) = M(6,7) = 

MC5,11) 

MC5,6) 

=s 

-^2^2 ^2b 



mC5,8) 

= 

Ejlj Sln2r 



MC5,9) 


-®1 =1 ■'19 SlnKj^g 



MC5,10) 

= 

*^2^2 ^20 



MC6,ll 


MC6,2) = M(6,3) = M(6,4) 

= mC6,9) = M(6,10) 

= M(6,ll) 

MC6,55 


-E^I^ SlnK^^ Ij 



mC6,61 

=: 

-^Zh 4b °°^hbh 



mC6,7) 


^zh ■'ib Si"** ^sb ■■2 



MC6,8) 

= 

^2h 4b ''2b h 



MC7,ll 


MC7,2) = MC&,3) = mC7,4) 

= MC7,9) = MC7,10) 

= M(7.li:) 

MC7,5) 

=: 

-Vz4b ''2bl2 



MC7,6) 


Egl^ 4^ Sin Ij 



MCr,7) 

rs 

®2^2 ^2b ^2b ^2 



MC7,8) 

= 

4^ Sink I 2 
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=: 

MC8,1) = MC8,3) = M(8,4) = M(8,5) = M(8,6) = M(8 


:= 

M(8,8) = M(8,9) = 0 

?^C8,10) 

= 

^2 ^2 ^20 ^20 ^2 

mC8,11) 


-G^ Sin ±2 

MC9,1) 

=: 

Sin 1-j^ 

M(9,2) 

- 

Cos 1;^ 

MC9,3) 

= 

Slnh 

mC9,^5 

= 

Cosh 

m(9,55 

= 

MC9,7) = MC9,9) = MC9,10} = MC9,11) = 0 

MC9,6) 

=: 

-1 

MC9,8) 

r: 

~1 


=: 

^Ib ^1 

M(10,2) 


Sin 1, 

MC10,31 

= 

Kib Cosh l 3 _ 

MC10,4) 

s 

Sinh 1^ 

m(io,5) 

= 

Cos2r 

M(10,6) 


MC10,8) = MC10,9) = MC10,10) = 0 

M(10,7) 

= 

Cos 2r 



M(ll,2) = MCll,3) = MCll,4) = MC11,6) = MC11,8) 



Mai,10) = 0 


= 

Sin2r 

mC 11,-7). 


Sln2r 

MC11>9) 

=: 

Cos 1;^ 


Cos2r 





4 


Date Slip 

Ycf23^hlj book Is to be returned on the 
date last stamped. 



fy]E- nsb-M- 



